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ABSTRACT 



To date, mid-infrared properties of Galactic black hole binaries have barely 
been investigated in the framework of mult i- wavelength campaigns. Yet, stud- 
ies in this spectral domain are crucial to get complementary information on the 
presence of dust and/or on the physical processes such as dust heating and ther- 
mal bremsstrahlung. Here, we report a long-term multi-wavelength study of the 
microquasar GRS 1915+105. On the one hand, we aimed at understanding the 
origins of the mid-infrared emission, and on the other hand, at searching for cor- 
relation with the high-energy and/or radio activities. We observed the source 
at several epochs between 2004 and 2006 with the photometer IRAC and spec- 
trometer IRS, both mounted on the Spitzer Space Telescope. When available, 
we completed our set of data with quasi-simultaneous RXTE / INTEGRAL high- 
energy and/or Ryle radio observations from public archives. We then studied 
the mid-infrared environment and activities of GRS 1915+105 through spectral 
analysis and broad band fitting of its radio to X-ray spectral energy distribu- 
tions. We detected polycyclic aromatic hydrocarbon molecules in all but one 
IRS spectra of GRS 1915+105 which unambiguously proves the presence of a 
dust component, likely photoionised by the high-energy emission. We also argue 
that this dust is distributed in a disc-like structure heated by the companion 
star, as observed in some Herbig Ae/Be and isolated cool giant stars. Moreover, 
we show that some of the soft X-ray emission emanating from the inner regions 
of the accretion disc is reprocessed and thermalised in the outer part. This leads 
to a mid-infrared excess that is very likely correlated to the soft X-ray emission. 
We exclude thermal bremsstrahlung as contributing significantly in this spectral 
domain. 
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1. Introduction 

In a multi-wavelength study of microquasars, the infrared presents a particular interest 
as the accretion disc, the jets, or the companion star may all be detected. Nevertheless, 
most of the previous studies focused on the accretion— ejection phenomena as seen is 
X-ray /near- infrared/radio correlations, and did not take the mid-infrared (MIR) emission 
into account (see e.g. Mirabel et al. 1998; Ueda et al. 2002; Corbel & Fender 2002; Corbel 
et al. 2003; Chaty et al. 2003; Homan et al. 2005; Chaty & Bessolaz 2006; Russell et al. 
2006). Yet, getting both MIR photometric and spectroscopic information is crucial to 
investigate the presence of dust, the disc illumination, thermal bremsstrahlung from the 
accretion disc's wind, or the contribution of relativistic ejecta. 

1.1. GRS 1915+105 

Discovered by the WATCH all-sky X-ray monitor on board the GRAN AT satellite, on 
1992 August 15 (Castro-Tirado et al. 1992, 1994), GRS 1915+105 is the first microquasar 
in which apparent superluminal radio ejecta were detected (Mirabel & Rodriguez 1994). 
The nature of its companion star was the subject of debate until Greiner et al. (2001b) 
unambiguously showed that it was a K/M red giant by detecting CO absorption features in 
its near-infrared (NIR) spectrum. Moreover, the orbital period of the system and the mass 
of the compact object were found to be 33.5 ±1.5 days (recently refined to 30.8 ± 0.2 day, 
Neil et al. 2007) and 14 + 4 M , respectively (Greiner et al. 2001a). The inclination is 
66° ± 2°, and estimates of the distance fall in the range 6—12 kpc (Chaty et al. 1996; Fender 
et al. 1999; Chapuis & Corbel 2004). 

GRS 1915+105 is strongly variable, on time scales from seconds to days. Using 
extensive RXTE timing observations, Belloni et al. (2000) showed that its X-ray behaviour 
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could be divided in 12 distinct luminosity classes (up to 14 today, Klein- Wolt et al. 2002; 
Hannikainen et al. 2003), and that the source was carrying out transitions between three 
canonical spectral states, labelled A, B (strong disc domination), and C (corona-dominated, 
no disc). Previous multi- wavelength studies showed the existence of a strong connection 
between the accretion disc instabilities and plasma outflows. In particular, discrete ejecta, 
emitting through optically thin synchrotron, are believed to be triggered during the 
transition between the C and A states, expanding adiabatically in the environment and 
detectable gradually from the NIR to the radio domains (see e.g. Mirabel et al. 1996; 
Fender et al. 1997; Mirabel et al. 1998; Eikenberry et al. 1998, 2000; Rodriguez et al. 
2008a,b). Moreover, the other known class of radio ejecta — the compact jets, emitting 
simultaneously from the radio to the NIR through optically thick synchrotron — are only 
detected in the x luminosity class, which is only seen in the C state. The presence of such 
a jet is characterised by a flat spectrum with a roughly constant flux density between 50 
and 100 mJy Long periods of the x luminosity class during which compact jets are present 
are called plateau, and often precede or follow a giant ejection (see e.g. Foster et al. 1996; 
Pooley & Fender 1997; Fender et al. 1999; Dhawan et al. 2000; Klein- Wolt et al. 2002; Fuchs 
et al. 2003b). 



1.2. Previous MIR studies of microquasars 

Koch-Miramond et al. (2002) presented an ISO spectrophotometric study of 
Cygnus X— 3 in quiescence, and found the MIR continuum to be due to free-free emission 
from the winds of the Wolf-Rayet companion star with perhaps a contribution from a cold 
dust component (see Fuchs et al. 2006, for similar conclusions on SS 433). Moreover, Fuchs 
et al. (2003a) reported ISOCAM photometric data of the GRS 1915+105 obtained at two 
different epochs, during a flaring activity and a plateau state. On the one hand, they showed 
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that, despite strong uncertainties, the MIR flux of GRS 1915+105 had likely increased 
between the two observations, and they argued that during the plateau, the MIR emission of 
the source was likely due to a compact jet, without excluding bremsstrahlung. They, on the 
other hand, excluded the contribution of dust, which is the opposite of the conclusion reached 
by Muno & Mauerhan (2006) to explain the MIR excess they detected in the Spitzer/IRAC 
SEDs of A0620-00 and XTE Jl 118+480 in quiescence. Finally, Migliari et al. (2007) argued 
that X-ray/UV irradiation of the disc in the thermal state, and compact jet in the hard state 
might be responsible for the excess they detected in the MIR emission of GRO J1655— 40; the 
same conclusion was reached by Gallo et al. (2007) concerning the quiescence of A0620— 00, 
V404 Cyg, and XTE J1118+480, claiming that dust component is not statistically necessary. 

In this paper, we report a long-term multi-wavelength study of GRS 1915+105 focusing 
on its spectroscopic and photometric MIR emission. It aimed at understanding its origins, 
as well as its possible connection with the high-energy and radio domains. The observations 
and the data analysis are presented in Sect. 2, while Sect. 3 and Sect. 4 are devoted to the 
analysis of the broad band X-ray to MIR SEDs of the source built with high-energy and 
Spitzer data. We discuss the outcomes in Sect. 5 and we give our conclusions in Sect. 6. 

2. Observations and data analysis 

We performed, between 2004 October 2 and 2006 June 5, 16 photometric and 
spectroscopic observations of GRS 1915+105 with the IRAC photometer and the IRS 
spectrometer (PI Y. Fuchs), both mounted on the Spitzer Space Telescope. Moreover, we 
completed our set of MIR observations with quasi-simultaneous high-energy observations, 
including INTEGRAL data already presented in great detail in Rodriguez et al. (2008a,b) 
(revolution #373), and several monitoring observations with RXTE (PI Morgan, the data 
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are immediately public). We finally also made use of observations of GRS 1915+105 obtained 
at 15 GHz with the Ryle Radio Telescope (PI G. G. Pooley). Table 1 lists all the data 
we used in this study. In the following, for each quasi- simultaneous X-ray/MIR/radio 
observation, we will refer to the integer part of the MIR observation date. 

2.1. Spitzer's MIR observations 

GRS 1915+105 was observed with IRAC at 3.59, 4.50, 5.80, and 8.00 pim. Each image 
was the combination of 24 sub-exposures of 10 s each, giving a total integration time of 
240 s in each filter. We performed photometry on the Basic Calibration Data (BCD) using 
the software MOsaicker and Point source Extractor (MOPEX) vl8.2.2. BCD data are raw 
data on which the Spitzer pipeline performs dark subtraction, multiplexer bleed correction, 
detector linearisation, flat fielding, cosmic ray detection, and flux calibration. We used 
MOPEX for pointing refinement, mosaicking, coaddition, and fluxes measurement through 
PSF fitting in a 3-pixel radius aperture. They were then scaled to a 10-pixel aperture 
using the aperture-correction factors as given in the IRAC manual 1 . The source was always 
detected, in all filters, and the absorbed fluxes are listed in Table 2. The uncertainties 
include the 3% systematic errors due to flux calibration instabilities (Reach et al. 2005). 

Spectroscopy was performed with IRS using the SL2 (5.20 — 7.70 yum), SL1 
(7.40 - 14.50 /im), LL2 (14.00 - 21.30 /im) and LL1 (19.50 - 38.00 /im) modules with 
the IRS Peak-up option for a better pointing accuracy. Total exposures times were set to 
120 s, divided in 2 sub-exposures in SL1 and SL2, and 300 s, divided in ten sub-exposures 
in LL1 and LL2. BCD data were reduced following the standard procedure given in the 



1 http://ssc. spitzer. caltech.edu/irac/dh/iracdatahandbook3. 0.pdf 
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IRS Data Handbook 2 . The basic steps were bad pixel correction, sky subtraction, as well 
as extraction and calibration (wavelength and flux) of the spectra — with the Spitzer IRS 
Custom Extraction software Spice v2.1.2 — for each nod. Spectra were then nod-averaged 
to improve the signal-to-noise ratio (SNR). In SL1 and SL2, GRS 1915+105 was always 
detected with SNRs good enough to allow spectroscopic features identification (SNR>10), 
as shown in Figure 1. Unfortunately, in LL1 and LL2, we never managed to detect the 
source. 



2.2. RXTE observations 

The RXTE data were reduced with the LHEASOFT v. 6.7. All data products were 
extracted from user's good times intervals (GTI). GTIs corresponded to times when the 
satellite elevation was greater than 10° above the Earth limb, the offset pointing less than 
0.02°, and proportional counter unit #2 was active. In order to identify the classes, we 
extracted Is resolution light curves in the 2-60 keV range and in the three energy bands 
defined in Belloni et al. (2000), from the Proportional Counter Array (PCA). These bands 
are 2.0-5.7 keV (channels 0-13, PCA epoch 5), 5.7-14.8 keV (channels 14-35), and 
>14.8 keV (channels 36-255). The colours were defined as HR1=5. 7-14.8/2.0-5.7 keV 
and HR2=14.8-60.0/2.0-5.7 keV. The shift of gain between the different epochs of PCA 
leads to different absolute values of the count rates, hardness ratios (HR) and position in 
the colour-colour (CC) diagrams, but the general shape of a given class is easily comparable 
to those of epoch 3 (Belloni et al. 2000), and therefore allowed us to easily identify the 
variability class in each observation (see Rodriguez et al. 2008a). 16 s resolution light curves 
were extracted from standard 2 data between 2.0 and 18.0 keV. All these light curves were 



2 http: / / ssc.spitzer.caltech.edu/irs/dh/dh32.pdf 
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corrected for background, using the latest PCA background models available for bright 
sources. Spectra, including those showing clear spectral variations, were then extracted and 
averaged over the entire observations for a better spectral fitting. 

2.3. INTEGRAL observations 

Our study makes use of the INTEGRAL Soft Gamma-Ray Imager (ISGRI, Lebrun 
et al. 2003) - the low energy detector of the Imager On-Board INTEGRAL (IBIS) - to 
cover the 18.0 to ~ 300.0 keV energy range, and the X-ray monitors JEM— X (Lund et al. 
2003) to cover the 3.0—30.0 keV one. The data reduction process can be found in Rodriguez 
et al. (2008a,b). 

2.4. MIR data dereddening 

The MIR data were dereddened from an optical absorption Ay = 19.6 ±1.7 (Chapuis 
& Corbel 2004) using the extinction laws given in Chiar & Tielens (2006). In their paper, 
the authors derived the A\/Ak ratio rather than the usual A\/Ay ones. To express the 
absorption ratio in a standard way, we assigned to Ak the value derived using the extinction 
law given in Fitzpatrick (1999) for the diffuse ISM (R y = 3.1), which is A K = 0.111 x A y . 
Although it has recently been demonstrated that the long-believed universal NIR extinction 
curve showed sharp variations depending on the line-of-sight (Fitzpatrick & Massa 2009, 
and references therein), these variations are rather observed outside of the Galactic plane 
(Galactic centre excluded), and our choice to fix Ry to the value of the diffuse ISM has no 
strong incidences in the dereddening of the GRS 1915+105 data. 

Moreover, if the authors give a universal expression up to 8.00 /im, they then propose 
two laws that take silicate absorption at 9.70 and 18.00 /im into account. Both differ in the 
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sense that one is valid for the diffuse ISM and the other one for the Galactic centre. Our 
absorbed MIR spectra exhibit a strong silicate absorption feature at 9.70 /im and we tried 
to deredden each of them using both laws. The best results were systematically obtained 
with the modelling of interstellar extinction due to the diffuse ISM, the silicate feature in 
the Galactic centre being too strong. 

3. Results 
3.1. Light curves 

The fluxes listed in Table 2 show that GRS 1915+105 is strongly variable in all the 
IRAC filters, which confirms the variation detected by Fuchs et al. (2003a). The first step 
of our study was to compare the evolution of the source's MIR emission with its X-ray and 
radio ones. Fig. 2 displays the 3.59 /im (IRAC), 1.2-12 keV (ASM) and 15 GHz (Ryle) 
light curves of GRS 1915+105 covering the time interval from MJD 53200 to MJD 53900. 
Although we cannot claim for a correlation in the existing data, it is noticeable that 
the source is at its minimum level of MIR emission (or rather the minimum level of our 
measurements) when the microquasar has a steady and relatively low X-ray activity (around 
MJD 53280, MJD 53308, MJD 53857, MJD 53890). On the other hand, the maximum level 
of its MIR emission was measured when the source entered in a quite high and unstable 
X-ray activity, around MJD 53500 and MJD 53676. This behaviour could be a hint for an 
X-ray/MIR correlation. 

Moreover, the comparison with the Ryle telescope light curve could point towards a 
bimodal relation between the radio and the MIR activities. Indeed, each time GRS 1915+105 
is detected at 15 GHz with a flux higher than 22 mJy (MJD 53496, 53500, and 53676), the 
MIR emission appears to be at its highest, while it is at its lowest when the radio flux is less 
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than 1 mJy (see Fig. 2 and Table 1). Nevertheless, it is important to point out that during 
all our spectroscopic and photometric observations, GRS 1915+105 has never been found 
in the plateau state; a MIR emission from the compact jets can consequently be excluded 
from now on. Each time the Ryle flux is high, the source is rather in the decaying phase of 
a giant discrete ejection, whose emission is optically thin. This MIR/radio bimodal relation 
is therefore very likely related to the disc activity. 



3.2. RXTE I INTEGRAL SEDs 

The quasi-simultaneous high-energy spectra were fitted in ISIS v. 1.5.0, between 3.0 
(6.0) and 20.0 keV for PCA (JEM-X), as well as 20.0 and 200.0 keV for ISGRI and 
HEXTE. In most cases, they were further rebinned so as to obtain good quality spectra. 
We added 3 % systematic errors to the uncertainties in the JEM-X spectra, 2% to the 
ISGRI ones and 1% to the PCA ones before the fitting process. In all fits a normalisation 
constant was added to account for uncertainties in the cross calibration of the instruments. 

The model we used during the fitting processes consisted in the combination of a 
multicolour black body (diskbb, Mitsuda et al. 1984), accounting for the accretion disc 
emission, a comptonisation component (comptt, Titarchuk 1994) for the corona and a 
Gaussian for the iron feature at 6.4 keV when present, all modified by photo-electric 
absorption (phabs). The column density was fixed at 3.5xl0 22 atoms cm -2 , as measured 
by Chapuis & Corbel (2004) from radio observations. If this value is consistent with 
some measurements derived from X-ray observations (see e.g. Ebisawa 1998; McClintock 
et al. 2006), several authors derived higher column densities clustered in the range 
5-8xl0 22 atoms cm" 2 (see e.g. Klein- Wolt et al. 2002; Lee et al. 2002). Nevertheless, our 
attempts with higher values always gave equally good fits and did not change significantly 
the parameters of the disc and the comptonisation components. Our goal being to fit broad 
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band SEDs including both X-ray and MIR data, we decided to use the same column density 
value in both spectral domains. 

A discussion on the validity of the use of comptt to model the hard X-ray emission can 
be found in Rodriguez et al. (2008b). Following the authors, the optical depth was fixed to 
0.01, the lowest allowed value, when GRS 1915+105 was in the states A or B, during which 
the hard X-ray emission is better fitted by a power law. We made this choice of comptt 
because the use of a power law strongly overestimates the contribution of the corona at 
low energy and consequently underestimates the disc's one, which has consequences in the 
optical and infrared domains. Moreover, a power law diverges once the data are dereddened, 
which forbids a spectral fitting from the X-rays to the MIR. comptt on the contrary peaks 
at 3kT disc and is negligible in the optical and the infrared. 

Finally, the Gaussian width was fixed at 0.8 keV and the iron line feature's energy was 
allowed to vary between 5.0 and 7.0 keV. Table 3 gives the best-fit parameters found for 
the spectra obtained on MJD 53280, MJD 53284, MJD 53500, MJD 53676, MJD 53851, 
and MJD 53890: kT disk and Norm are the accretion disc's temperature and norm, and kT e 
and r are the hard component's electrons temperature and opacity. The spectra obtained 
on MJD 63660 and MJD 53661 did not need any disc component but an extra power law. 
Their best-fit parameters are listed in Table 4. 



3.3. X-ray to MIR SEDs: the photometric case 

The first step to model the MIR emission of GRS 1915+105 was to understand 
in which extent the IRAC fluxes could be explained by the Rayleigh- Jeans tail of the 
accretion disc combined with the stellar emission. We then built the X-ray to MIR SEDs of 
GRS 1915+105 with the RXTE / INTEGRAL and IRAC data obtained quasi-simultaneously 
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on MJD 53284 and MJD 53890, i.e. when the measured photometric fluxes of the source 
were at their lowest. Concerning the latter date, it is worth noting that about 1.2 days 
separate the IRAC data from the RXTE ones. Nevertheless, the light curve displayed in 
Fig. 2 shows that between MJD 53890 and MJD 53892, the source had a steadily increasing 
1.2—12.0 keV X-ray emission (with no flares) and that the disc parameters were therefore 
barely variable. 

For each SED, the dereddened fluxes were stored in an ASCII file. They were then 
fitted along with the quasi- simultaneous high-energy data into ISIS, combining a spherical 
black body (bbodyrad) and the model described above. Greiner et al. (2001b) showed 
that the companion star of GRS 1915+105 was a K/M giant whose temperature is about 
T* = 4800±5oo K, and Harlaftis & Greiner (2004) found that its radius was around 19 R Q . 
During the fitting process, we then considered two distinct cases: 

• first, the temperature was allowed to vary between 2800 to 5000 K, which is the 
temperature scale of K/M giant stars (see e.g. van Belle et al. 1999). We obtained 
good fits, but for temperatures that systematically pegged at the minimum allowed, 
and for stellar radii clustered in the range 38—77 R Q depending on the considered 
distance. This is very unlikely because a giant star with such a low temperature would 
rather have a radius > 150 R Q (Dumm & Schild 1998; van Belle et al. 1999), 

• second, the star's temperature and radius were fixed to those of a K2 giant star, i.e. 
4520 K and 21 R Q as given in van Belle et al. (1999), and the source's distance was 
allowed to vary between 6 and 12 kpc. In both cases, the best fits, obtained for a 
distance of 6 kpc, were unable to completely reproduce the IRAC fluxes as there 
always was a MIR excess (see the fits displayed in Fig. 3). 

We then conclude that even when it is at its lowest, and in absence of any radio activity, 
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the MIR emission of GRS 1915+105 cannot be explained only by the companion star and 
Rayleigh- Jeans tail of the accretion disc. Moreover, on MJD 53500 and MJD 53676, the 
IRAC fluxes of the source almost doubled compared to their lowest values. It is impossible 
that this increase comes from the companion star. On the contrary, the disc dominates 
the 3.0—200.0 keV unabsorbed X-ray emission, which was not the case on MJD 53284 and 
MJD 53890 (Table 3). This therefore suggests a relation between the MIR and the disc 
activities of GRS 1915+105. 

4. Origin of the MIR excess of GRS 1915+105 

Although we excluded contributions from the companion star and the compact jets, 
the MIR excess we detected as well as the increase of the MIR fluxes of GRS 1915+105 may 
still be explained by (1) the presence of a dust component, maybe heated by the companion 
star and/or the X-ray/UV emission, (2) thermal bremsstrahlung from the accretion disc's 
winds, (3) illumination of the accretion disc, (4) optically thin synchrotron from a discrete 
ejection. 

4.1. A photoionised dust component? 

All the IRS SL1/2 absorbed spectra of GRS 1915+105 are displayed in Fig. 1, and 
Table 5 lists all the features we detected. All the measurements were carried out using 
the task IDEA of the data reduction software SMART v. 6.4.0. The central wavelength Agt, 
equivalent width W, full-width at half-maximum FWHM, flux, and signal-to-noise ratio 
SNR of each feature were computed through gaussian fitting. The highest source of error 
on the equivalent width and line flux measurement is due to the continuum, particularly 
uncertain at these wavelengths. Each time, instead of trying to fit it globally, the continuum 
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was assessed in the vicinity of each feature through a linear function fitting. Several 
attempts showed that the resulting systematic errors were about 5% of the measured flux, 
which were quadratically added to the statistical uncertainties. 

Along with several H I and H II emission lines that likely originate from the accretion 
disc, and a strong silicate absorption feature, we detected in each spectrum but the one 
obtained on MJD 53511, the so-called unidentified infrared features at 7.70 /im (with 2 
primary components at 7.60 and 7.80 /xm) and 11.25 /im (detected between 11.20 and 
11.30 //m). These lines are thought to be created by the family of the polycyclic aromatic 
hydrocarbon molecules (PAH, Leger & Puget 1984; Puget et al. 1985; Allamandola et al. 
1985), which are found in the MIR spectra of many objects with associated dust and gas 
components illuminated by UV photons. 

The presence of such features in the MIR spectrum of GRS 1915+105 infers that 
there is dust in the system. Moreover, despite the strong uncertainties, all the PAH lines 
appear to be strongly variable in flux; this might be related to the photoionisation state 
of the environment of GRS 1915+105. Indeed, the excitation mode of PAH molecules 
is vibrational; the 7.70 /iin feature is due to C— C stretching, whereas the 11.25 /im is 
due to C— H out-of-plane bending. The C— H bond being weaker than the C— C one, the 
photoionised PAH is principally detected at 7.70 /im. The flux ratio F7.7/F11.3 is therefore a 
good indicator of the photoionisation state of the environment in which PAH molecules are 
detected, and the larger it is, the higher the photoionisation degree (see e.g. Allamandola 
et al. 1989; Tielens 2008, for comprehensive reviews on PAH properties). 

The companion star of GRS 1915+105 is a late K giant. If its emission can heat up 
a potential dust component, it is impossible that it photoionises it as it is too cold to 
emit enough energetic UV photons. Therefore, only the X-ray emission can be responsible 
for the photoionisation of the PAH molecules. To validate this hypothesis, we compared, 
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for each spectrum, the evolution of the F 7 , 7 /Fu_ 3 PAH flux ratio with the corresponding 
1.2-12.0 keV X-ray flux of GRS 1915+105. The result is displayed in Fig. 4. Despite the 
strong uncertainties due to the silicate absorption, it suggests a correlation; the dust may 
therefore be photoionised by the X-ray/UV photons originated from the accretion disc 
and/or the corona. This would partly explain why no PAH molecules are detected in the 
GRS 1915+105 MIR spectrum obtained on MJD 53511. Indeed, a [Ne II] emission feature 
is present at 12.81 /im. The Ne atom ionisation potential is about 21.56 eV, while the PAH 
one is smaller than 10 eV, depending on the molecule size (Ruiterkamp et al. 2005). The 
absence of the 11.25 /im PAH feature could consequently mean that almost all the PAH 
was photoionised. But this spectrum also displays a huge increase of the continuum below 
about 9.00 /im, reaching the flux level measured through photometry on MJD 53500 i.e. 
about 8.00 mJy and 5.00 mJy at 5.80 /im and 8.00 /im, respectively. The reason for this 
increase will be discussed in the next section, but the non-detection of the PAH feature at 
7.70 /im could be due to a contamination by the continuum, as the peak flux of the feature 
appears to be less than 6.00 mJy in all the other spectra where it is detected. 



4.2. Effect of the thermal bremsstrahlung 

Thermal bremsstrahlung from an expanding wind might be partially responsible for the 
MIR emission of GRS 1915+105. Indeed, Begelman et al. (1983) showed that X-ray driven 
winds could form above an accretion disc heated by X-ray radiation with luminosity a few 
percents above the Eddington limit, which is likely the case of GRS 1915+105 as Lee et al. 
(2002) and more recently Neilsen & Lee (2009) and Ueda et al. (2009) detected such wind 
in the system. 

Thermal bremsstrahlung was already invoked in van Paradijs et al. (1994) to explain 
the strong MIR excess the authors had detected in the emission of GRO J0422+32. In their 
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paper, they used the formalism given in Rybicki & Lightman (1979) to assess the expected 
10.80 /iin luminosity of a spherical expanding wind. We followed the same steps, for a 
disc wind emitted in a solid angle Q, to assess the expected monochromatic luminosity of 
GRS 1915+105 due to free-free emission at 8.00 pim. In such a wind, the mass-loss rate M w 
is: 

^ w = ( 7~ J x 47rr 2 m p n e w w (1) 



where r is the distance within the wind, m p the proton mass, n e the electronic density, 
and v w the wind velocity. Following Rybicki & Lightman (1979), the thermal bremsstrahlung 
emissivity at the frequency v can be written as: 

«r ( Ti, 2 \ __hv 

t v = 6.80 x 10~ 45 l-^= J e W x g W cm" 3 Hz" 1 (2) 

where T is the wind temperature, g the Gaunt factor, and h and k-Q the Planck and 
Boltzmann constants, respectively. In the following, we fix g to 1, corresponding to a large 
angle regime in the interaction between an electron and an ion. 

Replacing n e from Eq. 1 into Eq. 2, and integrating e v over the radial distance r 
between the launching radius R (in cm) and infinity, we obtained the monochromatic 
luminosity L u at the frequency v. 

L u = (J^J J e u x 4nr 2 dr (3) 
= 2.04 x 10 s e —= . ' " WHz" 1 (4) 

Following Ueda et al. (2009), we can approximate the mass-loss rate M w as: 

M w w 1.00 x 10 11 (J^J v w kg s" 1 (5) 



which in turn gives the following expression for the monochromatic luminosity: 

hv 

L v = 2.04 x 10 30 ( ^= (6) 



V4tt/ R Vr 
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In their thermally-driven wind model, Begelman et al. (1983) introduced four important 

parameters, the Compton temperature T c , where heating from Compton scattering and 

cooling from inverse Compton are balanced out, the Compton radius R c , for which the 

escape velocity equals the isothermal sound speed at the Compton temperature, the critical 

luminosity L CT , above which a Compton heating disc wind can overcome gravity, and T c h, 

which is the characteristic temperature of a parcel of gas that rises at a height Rq above the 

disc in a finite heating time. Those parameters are defined as: 

1 f UN hvL v dv 

Tc = IF %r 1 K ( 7 ) 
4fc B J Ui L v dv 

9.80xl0 17 M x 

^ C = T C Me Cm (8) 
L cr « 2.88 x 10 2 -p= (9) 

2 _ 2 

L V / R \ s 



r * = rc l^J (10) 

where M x is the black hole's mass, L the X-ray bolometric luminosity, and L E the 
Eddington luminosity. On MJD 53284, using the continuum parameters given in Table 3, 
and integrating the system's monochromatic luminosity between 1.0 and 1000.0 keV, we 
find T c w 5.80 x 10 6 K, which leads to R c « 2.37 x 10 12 cm, T cr w 0.12 x L E , and 
L w 0.37 x L E . 

Begelman et al. (1983) and Woods et al. (1996) showed that a disc wind could 
develop for R > 0.2 x R c and would get strong for R > R c . It is reasonable to assert 
that the thermal bremsstrahlung arises from the gravity-free part of the wind, in which 
T Rj T ch < T c . Then, considering L w 3 x T cr , Eq. 10 leads to Rq > 3 x R c . Replacing 
T by Tc and Rq by 3 x i?c in Eq. 6, the inferred luminosity at 8.00 /im, L$ is: 

L 8 w 1.19 x 10 14 W Hz" 1 (11) 

On MJD 53284, the unabsorbed luminosity of GRS 1915+105 at 8.00 /im is about 
8.27 ± 1.30 x 10 13 W Hz -1 . With a covering factor of about 5%, as measured for 
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GRS 1915+105 in Neilsen & Lee (2009), we deduce a bremsstrahlung -induced luminosity at 
8.00 jum of about (5.60 ± 0.55) x 10 12 W Hz -1 , which is an order of magnitude lower that 
what we measured. Even for a larger value of the covering factor as high as 0.2 (Proga et al. 
2000), we derive about (2.38 ± 0.23) x 10 13 W Hz" 1 , which is still too low. It is therefore 
very likely that thermal bremsstrahlung from the accretion disc wind barely contributes 
to the MIR flux of GRS 1915+105. Moreover, Fig. 1 shows that on MJD 53511, the MIR 
continuum of the source below about 9 /zm increased while it remained almost the same 
beyond. Yet, the spectral signature of free- free emission is a power law with an index from 
(optically thin) to 2 (optically thick, Wright & Barlow 1975), the wind being in the optically 
thin regime in the infrared. An increase of the MIR flux due to bremsstrahlung would be 
detected at all wavelengths, which then excludes it as a reason for the MIR brightening. 

4.3. X-ray to MIR SEDs: irradiation of the disc 

In the soft state, reprocessing — in the outer disc — of X-ray and UV photons originated 
from the inner part likely dominates the UV and optical emission of microquasars (Vrtilek 
et al. 1990; Fukue 1992; Sanbuichi et al. 1993; van Paradijs & McClintock 1994; Hynes 
et al. 1998; Esin et al. 2000; Hynes et al. 2002). In the hard state, there might also be a 
contribution of the reprocessed hard X-ray photons from the corona. Indeed, Ueda et al. 
(2002) showed that it could even represent about 20—30 % of the K-b&nd emission of 
GRS 1915+105 in the plateau state. 

We showed that the sharp increase of the MIR emission on MJD 53511 — only detected 
below about 9.00 /im- was not due to thermal bremsstrahlung. Dust heating is also 
little plausible as the MIR increase would have been detected at all wavelengths. On the 
contrary, it seems rather consistent with irradiation of the outer disc and/or optically thin 
synchrotron from a discrete ejection. To confirm this hypothesis, we built — for MJD 53284, 
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MJD 53511, and MJD 53851 - the X-ray to radio SEDs of GRS 1915+105, including the 
dereddened MIR spectra. Moreover, on MJD 53284 and MJD 53851, GRS 1915+105 was 
not detected in the radio domain so we did not include any radio flux in these SEDs. 

The way we built the MJD 53511 radio to X-ray SED of GRS 1915+105 deserves some 
justifications. In the high energy domain, we used the MJD 53500 RXTE data because 
we did not have any high-energy observations quasi- simultaneous with the IRS ones (see 
Table 1). Indeed, the ASM fluxes and C/A hardness ratios being similar at both dates 
(see Fig. 5), it is likely that the disc parameters were the same. Moreover, to compensate 
the lack of data in the radio domain, we made use of the archival VLA fluxes of the 
source obtained on MJD 53513 at 8.46, 14.94, and 22.46 GHz (37.90, 23.20, and 16.20 mJy 
respectively, about 10% uncertainties) 3 . At this epoch, GRS 1915+105 was in the decaying 
phase of a giant ejection, and the VLA flux at 14.94 GHz is similar to the ones from the 
Ryle telescope on MJD 53510 and MJD 53512, proving that the unknown MJD 53511 radio 
fluxes are barely different from the VLA ones. 

The model we chose to take the reprocessing into account is diskir (Gierlihski et al. 
2008, 2009). Roughly, it is a extension of diskbb that includes disc irradiation (both from 
the inner region and the corona) as well as comptonisation (based on nthcomp). The model 
has nine parameters: the disc's temperature /cTdi SC and norm Norm (same as diskbb), the 
hard X-rays power law T and temperature kT e , the ratio between the corona's and the disc's 
luminosity L c / L d , the fraction of hard X-ray emission that illuminates the disc f in , the 
irradiated radius R irr expressed in term of the disc inner radius, the fraction of soft X-ray 
emission which is thermalised in the outer disc / ou t, and the logarithm of the outer radius 
Irout expressed in function of the inner radius. The first seven parameters are completely 
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defined by the high-energy data while the two latter are characterised by the optical and 
infrared ones. 

diskir was additively combined to a gaussian accounting for the iron feature at 6.4 keV 
(0.8 keV frozen width), and we fitted all the data, in ISIS, in three steps: 

• we first modified the model with a photo-electric absorption phabs (Ah = 

3.5 x 10 22 atoms cur 2 ) and fitted the high-energy data only. The electron temperature 
was frozen to the value found in the previous fits with comptt (see Table 3), the 
irradiated radius was frozen to 1.1 x R m after several unsuccessful attempts that 
showed that it was poorly constrained, f\ n was fixed to 0.1 for the corona-dominated 
spectrum and to 0.3 for the disc-dominated one (for a fixed 0.1 disc's albedo, Poutanen 
et al. 1997; Ibragimov et al. 2005; Gilfanov 2009), and / out and Irout were frozen to 
and 3, respectively, as the X-ray data do not allow to constrain them, 

• we built the new radio to X-ray SEDs with the unabsorbed MIR spectra — stored in 
ASCII files - the RXTE data sets and the VLA radio fluxes (MJD 53511 only), 

• we finally fitted the global SEDs with the previous model combined to spherical 
black body component accounting for the companion star emission, another one 
accounting for the detected dust emission, and a power law for the radio emission 
(MJD 53511 only). All the parameters were allowed to vary freely, except the stellar 
ones (temperature and radius fixed to 4520 K and 21 R Q , respectively), kT^c and i?; rr 
(same values as in step 1), and Irout which was frozen to 0.52a (Chaty et al. 2003), a 
being the orbital separation derived from the third Kepler's law for a 30.8 days orbital 
period (Neil et al. 2007), and a 14 M & and 0.86 M & black hole and companion star, 
respectively (Harlaftis & Greiner 2004). Note that Irout was dynamically tied to Ri n 
for a 10 kpc distance and a 66° inclination. 
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All the best-fit parameters are listed in Table 6, and the fitted SEDs are displayed in 
Fig. 6. The lack of optical and NIR fluxes being a caveat to constrain irradiation, we added 
in each SED the if-band magnitudes of GRS 1915+105 found in Neil et al. (2007), the 
uncertainties corresponding to the magnitudes the day before and after our observations. 
These data were not used to fit the SEDs, but they show the emission level given in the 
i^-band by the irradiation component is not inconsistent with the measurements at the 
same epoch. 



5. Discussion 

5.1. Origin of the dust component 

Our results suggest the existence of a cold dust component (T dust « 300 — 500 K) in the 
vicinity of GRS 1915+105 that likely interacts with the high-energy emission of the black 
hole binary. Moreover, the dust's temperature and radius appear to be roughly constant 
whatever the level of the X-ray emission, which implies that the dust is heated by the K 
giant companion star rather than the high-energy photons. And lastly, the average dust 
extension derived from the fits, -Rdust ~ 500 R Q or 3.5 x 10 13 cm, is roughly 5 times larger 
than the orbital separation and 10 times larger than the outer radius of the accretion disc, 
which means that the MIR emission due to the dust is produced well beyond the binary 
orbit and that the dust enshrouds the whole system. 

Many red giant stars are known to be embedded in a dusty shell that originates 
from the slow and dense stellar winds (see e.g. Hagen 1978; Zuckerman & Dyck 1986; 
Morris 1987; van Loon et al. 2005), and this might be the case of the companion star of 
GRS 1915+105. To check the consistency of the dust's temperature derived from the fits 
with the dusty stellar winds hypothesis, we can use the simple relation giving the expected 
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temperature of a spherical dust shell in a thermodynamic equilibrium with a central star. 
Following Rahoui et al. (2009), this temperature is: 



where h and k are the Planck and the Boltzmann constants, T the gamma and ( the 
Riemann zeta functions, and Q and n such as the chromatic grain emissivity Q u is defined 
as Q v = Qov n . Typical values for carboneous dust are n = 1.2 and Q = 1.52 x 10~ 8 r g , 
where r g is the average dust grain's radius (see e.g. Draine & Lee 1984; Robberto & Herbst 
1998). For 0.01 fim< r g <0.1 fim, as observed for interstellar dust grains (Draine & Lee 
1984), T* = 4520 K, = 21 R Q , and -Rdust ~ 500 R Q , the inferred dust temperatures is 
about 315 K < T dust < 490 K, which is consistent with our fits. 

However, the dust could also originate from a dusty disc-like circumstellar component. 
Such discs have already been invoked around some isolated first-ascent red giant stars, and 
a possible explanation for their presence could be the engulfment of an hypothetic low-mass 
companion when the star entered into the red giant phase (see e.g. Jura 2003; Jura et al. 
2006; Melis et al. 2009). They also have been detected around cataclysmic variables (Dubus 
et al. 2004; Howell et al. 2006; Brinkworth et al. 2007; Hoard et al. 2009), and Muno 
& Mauerhan (2006) suggested their presence around A0620-00 and XTE J1118+480 to 
explain the 8.00 /iin MIR excess in the emission of both sources while in quiescence. The 
two most common accepted explanations for the presence of circumbinary discs (CBDs) 
are (1) the dust was ejected from the binary with angular momentum during the common 
envelope phase, (2) it comes from a supernova fallback, as it was recently argued for the 
anomalous X-ray pulsar 4U 0142+61 (Wang et al. 2006). 

In the case of a flat and optically thick CBD irradiated by the companion star, the 




(12) 
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expected temperature at a radius r is (Chiang & Goldreich 1997): 

^'"-©'(t)' 5 (13) 

For a CBD whose inner parts are truncated by tidal forces, the expected minimum inner 
radius is about r in 1.7a, where a is the binary separation (Artymowicz & Lubow 1994). 
In the case of GRS 1915+105, this leads to a maximum inner radius temperature of an 
hypothetical flat and optically thick CBD of T CBD (r in ) ps 620 K, which again is consistent 
with the MIR excess due to dust that we detected. 

At this point, our data do not allow us to discriminate between a spherical and a disc 
geometry, but simple considerations confirm the presence of dust heated by the companion 
star and enshrouding the black hole binary. Moreover, it is worth noting that the MIR 
spectra of GRS 1915+105 are strongly similar in shape to those of some pre-main sequence 
or Herbig Ae/Be stars, silicate absorption and PAH features included (see e.g. van Boekel 
et al. 2004; Sloan et al. 2005; Boersma et al. 2009; Berne et al. 2009, for such spectra). 
Such stars are known to exhibit a MIR excess due to the presence of an equatorial dusty 
disc within which PAH molecules are photoionised by the UV emission of the central star, 
and we suggest that this similarity strengthens the disc scenario for the dust distribution 
around GRS 1915+105. 

Finally, the silicate absorption feature due to the diffuse interstellar medium is strongly 
correlated to the optical extinction as Ay = (18.5 + 2) x T9.7 (Draine 2003), where T9.7 is the 
optical depth of the silicate absorption at 9.70 /im. Following Chiar et al. (2007), we took 
T9.7 = — £n(F 9 .7/F cont i nuum ), where F9.7 and Fcontinuum are the source and continuum's fluxes 
at 9.70 /im, respectively. The optical depth was computed on the spectrum obtained on 
MJD 53851, as this is the one for which X-ray reprocessing has the smallest contribution in 
the the MIR. The continuum was fitted using the ranges 5.20—7.00 /im and 13.00—14.50 /xm 
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(in order to exclude the contribution of the silicate absorption feature) with a second 
order polynomial. This method is strongly uncertain, especially concerning the continuum 
fitting, and the result should therefore be considered with caution. We nevertheless infer 
Ay = 20.04 ± 4.13, which is consistent with the value found by Chapuis & Corbel (2004) 
for the optical extinction in the line of sight of GRS 1915+105. We therefore conclude that 
the silicate absorption feature in the MIR spectrum of GRS 1915+105 is likely due to the 
diffuse interstellar medium. 

5.2. Importance of irradiation 

The results of the fitting of the three SEDs displayed in Fig. 6 suggest that not only 
X-ray irradiation of the disc dominates the UV to NIR emission of GRS 1915+105, but 
that it also extends to the MIR, where it is overcome by the dust component between about 
6 to 10 /im, depending on the accretion disc's flux. In the thermal state (MJD 53511), 
the MIR continuum is even strongly dominated by the thermalisation in the outer region 
of the soft X-ray emission from the inner parts. These results might once again emphasise 
the peculiarity of GRS 1915+105. Indeed, only Migliari et al. (2007) previously proposed 
X-ray irradiation of the disc to explain the MIR emission of GRO 1655—40 in the thermal 
state, and van Paradijs et al. (1994) excluded it for GRO J0422+32, arguing that the 
accretion disc was not large enough to have outer regions sufficiently cold to emit at these 
wavelengths. This is not the case of GRS 1915+105 as the latter exhibits a large accretion 
disc, with an assessed outer radius of about 2.86 x 10 12 cm. A strong contribution of the 
X-ray irradiation extending to the infrared was then expected, and Ueda et al. (2002) 
already argued that it was responsible for about 20—30% of the K— band flux in the hard 
state, even in presence of compact jets. 
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Disc illumination therefore provides a consistent explanation for the variations of the 
MIR continuum of GRS 1915+105. A caveat nevertheless forbids definitive conclusions. 
Indeed, X-ray reprocessing is thought to be the dominant contribution to the UV and optical 
fluxes of X-ray binaries, and any model needs these information to be well constrained. 
Their lack in the set of data we used to fit the SEDs is then a strong limitation to 
our interpretation, as the process could have artificially increased the emission at these 
wavelengths to better fit the MIR data. So, even if the fluxes given in Neil et al. (2007) 
are consistent with the ones derived from our model, only the information on at least the 
quasi-simultaneous J— , H-, and K— bands magnitudes could strengthen our conclusions. 

6. Conclusion 

We presented a multi-wavelength study of GRS 1915+105 whose outcomes suggest 
that, in the absence of discrete or continuous ejecta, the MIR continuum of the source 
is mainly due to the X-ray irradiation of the accretion disc and to a photoionised dust 
component. This might have consequences on the interpretation of the MIR emission of 
microquasars in presence of compact jets. Indeed, dust might be ubiquitous around isolated 
compact objects and X-ray binaries, because of mass transfer during the common envelope 
phase or material from supernova fallback. If so, compact jets could contribute less than 
expected at infrared wavelengths, with perhaps a cutoff frequency in the millimeter domain. 
To confirm our results, it is therefore crucial to increase the sample of microquasars and 
systematically observe them through MIR spectroscopy, as this is the only way to obtain 
firm information on both their environment and their continuum. In particular, studying 
microquasars whose variation time scales are longer than the GRS 1915+105 ones, and 
which do not exhibit such rapid transitions between spectral states would strongly facilitate 
the mult i- wavelength observations and would allow to reach definitive conclusions. 
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Table 1. Summary of all the Spitzer's data of GRS 1915+105 we made use of in this study 



IRAC IRS RXTE INTEGRAL Ryle Ryle fluxes (mjy) Class 





53280.275 


53280.241- 


53280.348 


53280.636- 


-53280.746 


<1 


P 


53284.176 




53284.177- 


53284.236 


53283.673- 


-53283.706 


<1 


P 




53299.226 






53298.587- 


-53298.663 


<1 




53308.741 








53308.811- 


-53308.841 


<1 






53484.099 






53484.027- 


-53484.169 


76-91 




53496.287 








53496.209- 


-53496.230 


37-43 




53500.571 




53500.326- 


53500.430 


53500.202- 


-53500.358 


22-32 


S 




53511.697 














53636.571 


















53660.072 


53659.982- 


53660.079 


53659.577- 


-53659.906 


27-47 


<t> 




53661.182 


53661.706- 


53661.856 


53660.602- 


-53660.881 


18-55 


<t> 


53676.070 








53676.247-53677.488 53675.583- 


-63675.641 


64-75 


X, M 




53851.419 


53851.307- 


53851.319 


53851.102- 


-53851.320 


<3 


X 


53857.587 


















53874.818 














53890.827 




53892.047- 


53892.063 


53888.239- 


-53888.261 


<1 


P 



Note. — For each instrument, we give the day of observation (in MJD), and, when available, we give the day of quasi-simultaneous 
coverage with RXTE, INTEGRAL, and/or the Ryle telescope, as well as the Ryle flux level (15 GHz) in mjy. Moreover, when 
high-energy observations were available, we give the spectral class of the source as defined in Belloni et al. (2000) 
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Fig. 1. — IRS spectra of GRS 1915+105 from 5.20 to 14.50 /xm. All detected features are marked. 
For the MJD 53511 spectrum, the arrows mark the position of the undetected 7.7 and 11.25 /jm 
PAH features. 
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Fig. 2.— Ryle (15 GHz, top), ASM (1.2-12.0 kev, middle), and IRAC (3.59 fxm, bottom) light 
curves of GRS 1915+105 between MJD 53200 and MJD 53900. 
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Fig. 3.- GRS 1915+105 X-ray to MIR SEDs - fitted with the model 

diskbb+gaussian+comptt+bbodyrad - built with the MJD 53284 and MJD 53890 
RXTE /PCA+HEXTE and IRAC data. A MIR excess is detected. 
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Table 2. Absorbed fluxes (in mJy) of GRS 1915+105 in the four IRAC filters 



MJD 


3.59 jum 


4.50 


5.80 


8.00 fj,m 


53284 


5.70±0.17 


5.32±0.16 


4.80±0.15 


3.12±0.10 


53308 


5.26±0.16 


4.85±0.15 


4.32±0.13 


2.95±0.10 


53496 


8.69±0.26 


8.12±0.24 


7.61±0.23 


5.03±0.16 


53500 


10.28±0.31 


10.09±0.30 


8.80±0.27 


6.14±0.19 


53636 


6.20±0.19 


5.88±0.18 


5.33±0.16 


3.77±0.12 


53676 


10.70±0.32 


10.31±0.31 


9.15±0.28 


6.41±0.20 


53857 


4.82±0.15 


4.75±0.14 


4.23±0.13 


2.87±0.09 


53890 


5.16±0.16 


4.92±0.15 


4.36±0.13 


3.05±0.10 



Note. — Uncertainties are given at la and include 3% sys- 
tematic errors. 



Table 3. Best parameters obtained from the fit of the MJD 53280, MJD 53284, MJD 53500, MJD 53676, MJD 53851, and 

MJD 53890 high-energy spectra of GRS 1915+105 





Parameters 


MJD 53280 


MJD 53284 


ATT TA roCAA 

MJD 53500 


A T TTA tTO/""7/"' 

MJD 53676 


MJD 53851 


A T TTA roonn 

MJD 53890 




fcT disc (keV) 

Norm 
kT c (keV) 


-, 97 +0.02 
1-z ' -0.03 

594.0^ 
158.8±1;5 


765.9±^; 8 
155.2± 3 ;° 


1 - oo -0.02 

267.5±l?;g 
118.lt 8 ;* 


2 05+ ' 04 

z - uo -0.04 

90.99±i;fl 
158. 2± 7 '; 8 


n 01 +0.05 
U - 8i -0.05 

1 1 93 n +386 - 6 

113 °- u -264 .6 

07 o + 13.7 
z ' - J -6.6 


-, O o+0.07 
1 - zo -0.09 

396.8±lf;i 
131.7™ 




T 


0.01 (frozen) 


0.01 (frozen) 


0.01 (frozen) 


0.01 (frozen) 


Q5+ - 37 


0.01 (frozen) 


-Ftotai a 


(xlCT 8 erg cm~ 2 s _1 ) 


3.08 


3.49 


4.64 


4.59 


1.71 


2.13 


-^disc 


(xl0~ 8 erg cm~ 2 s _1 ) 


1.22 


1.19 


3.75 


3.74 


0.20 


0.85 




X 2 (d.o.f) 


0.93 (53) 


1.09 (53) 


1.29 (53) 


1.65 (49) 


1.31 (53) 


1.25 (53) 



Note. - The best-fit model is phabs(diskbb+gaussian+comptt), and the errorbars are given at the 90% confidence 
level. These spectra were built with RXTE /PC A+HEXTE data, except on MJD 53676 for which we used INTE- 
GRAL/ JEM-X+ISGRI data. 

a -Ftotai and -Fdi sc are the total and disc unabsorbed fluxes, extrapolated to 3.0—200.0 keV. 
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Table 4. Best parameters from the fit of the MJD 53660 and MJD 53661 
RXTE /PCA+REXTE spectra of GRS 1915+105 



Parameters MJD 53660 MJD 53661 



kT sccd (keV) 


n 5s +0 10 


n 7 9+0.05 
u -' z -0.04 


kT c (keV) 


6 17+ 1 - 84 

°- 1 '-0.94 


5-48lS:S 


T 


Z.UO_ 42 


2.0418:3* 


pa 


9 71 +0.30 


2.87±°;i° 


^totai b (xl0- 8 crg cm- 2 s" 1 ) 


3.08 


4.43 


X 2 (d.o.f) 


0.82 (53) 


0.74 (53) 



Note. The best-fit model is 

phabs(comptt+gaussian+powerlaw), and the errorbars 
are given at the 90% confidence level. 

a r is the power law photon index. 

b -Ftotai is the total unabsorbed flux, extrapolated to 
3.0-200.0 keV. 
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Table 5. List of all detected features for each IRS spectrum of GRS 1915+105 



Features 


A 


Afit 


W 


FWHM 


Flux 


SNR 




H 


W 


H 


(H 


xl0~ 21 W cm- 2 




MJD 53280 


HI 


7.508 


7.522±0.004 


-0.041 


0.113±0.005 


1.030±0.119 


11.29 


PAH 


7.800 


7.801±0.008 


-0.032 


0.228±0.009 


1.266±0.137 


15.93 


PAH 


11.250 


11.251±0.004 


-0.193 


0.178±0.004 


0.622±0.041 


18.04 


He II ? 


12.367 


12.353±0.003 


-0.090 


0.163±0.003 


0.424±0.026 


27.54 


HI ? 


12.587 


12.590±0.011 


-0.075 


0.239±0.014 


0.576±0.089 


14.65 


MJD 53299 


HI 


7.508 


7.519±0.003 


-0.058 


0.134±0.004 


1.546±0.169 


34.30 


PAH 


7.600 


7.621±0.002 


-0.034 


0.100±0.003 


1.432±0.078 


25.37 


PAH 


11.250 


11.206±0.003 


-0.227 


0.189±0.003 


0.599±0.039 


36.64 


MJD 53484 


HII 


7.460 


7.458±0.004 


-0.050 


0.122±0.007 


1.546±0.06 


68.28 


PAH 


7.700 


7.722±0.008 


-0.097 


0.240±0.009 


2.581±0.190 


16.55 


PAH 


11.250 


11.299±0.014 


-0.110 


0.143±0.012 


0.463±0.060 


8.85 


HI 


12.587 


12.599±0.007 


-0.140 


0.194±0.007 


0.670±0.076 


12.49 


MJD 53511 


HII 


7.460 


7.460±0.002 


-0.029 


0.087±0.002 


0.903±0.073 


17.47 


[Nc II] 


12.813 


12.803±0.001 


-0.080 


0.136±0.001 


0.427±0.011 


60.68 


MJD 53660 


PAH 


7.700 


7.719±0.021 


-0.029 


0.264±0.017 


1.293±0.138 


16.11 


PAH 


11.250 


11.279±0.011 


-0.077 


0.256±0.005 


0.692±0.035 


46.56 


MJD 53661 


PAH 


7.800 


7.812±0.005 


-0.051 


0.187±0.006 


1.389±0.111 


14.66 


PAH 


11.250 


11.311±0.005 


-0.100 


0.208±0.005 


0.644±0.050 


17.89 
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Table 5 — Continued 



Features A A fit W FWHM Flux SNR 

(fim) (/mi) (fim) (fim) xlO~ 21 Wcm~ 2 

HI 12.157 12.136±0.004 -0.094 0.183±0.006 0.625±0.055 17.69 

HI 12.611 12.606±0.011 -0.089 0.223±0.012 0.667±0.099 11.21 

MJD 53851 

PAH 7.800 7.773±0.006 -0.043 0.119±0.006 0.853±0.126 11.29 

PAH 11.250 11.269±0.009 -0.347 0.281±0.010 1.188±0.120 11.42 

MJD 53874 

PAH 7.800 7.831±0.010 -0.026 0.121±0.010 1.152±0.112 11.76 

PAH 11.250 11.247±0.001 -0.076 0.115±0.002 0.662±0.022 45.99 

HI 12.587 12.568±0.011 -0.101 0.253±0.028 0.670±0.128 11.02 



Note. — We give the name of the feature, its laboratory (A) and measured (Afit) wavelengths, 
its equivalent width W, its full-width at half-length (FWHM), its flux, and its signal-to-noise 
ratio (SNR). 
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Fig. 4. — Ratio of the PAH lines absorbed fluxes at 7.70 and 11.25 fim in function of the 1.2- 
12.0 keV ASM absorbed flux. Errorbars are given at 3a. 
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Fig. 5.— Top: day-averaged ASM light curve of GRS 1915+105 between MJD 53495 and 
MJD 53520. 

Bottom: C/A ASM hardness ratio. 
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Table 6. Best parameters obtained from the fits of the GRS 1915+105 MIR to X-ray 

dereddened SEDs 



Parameters 


MJD 53280 


MJD 53511 


MJD 53851 


fcT disc (keV) 


i no+0.02 
1 - JO -0.02 


1 qo+0.03 


u - ' °-0.06 


Norm 


520.08^32-59 


251.46±g;?2 


1756.80t^-3i 


Irout (Ri n ) 


6.01 


6.17 


5.74 


r 


2.74lg;g 


2 96+ 07 

z - yD -0.07 


r> i O+0.02 

z - lo -0.03 


Lc/Ld 


o-4fi±8:gi 


nfi7+ 002 

U - UD ' -0.001 


1 QO+0-06 
J--o»_0 .05 


fin 


0.1 (frozen) 


0.3 (frozen) 


0.1 (frozen) 


fout (icr 4 ) 


29.39tg 1 7 f 


eo+i05.58 
ooo.oo_ 8g 2 Q 




T dust (K) 


396.55t&£ 


306.66±f£^ 


461.16l 8 6 6 ; 5 3 ° 


-Rdust (Rq) 


538.401^ 


743.66±^1;^ 


07c 1 7+87.91 
0(0.1/ _ 72.91 


Radio 




u - oo -0.09 




X 2 (d.o.f) 


1.28 (187) 


1.12 (190) 


1.64 (187) 



Note. The best-fit model is 

phabsfdiskir+gaussian+bbodyrad+bbodyrad+powerlaw), and 
the errorbars are given at the 90% confidence level. The SEDs 
were built with the MJD 53280, MJD 53511, and MJD 53851 
RXTE /PCA+HEXTE and IRS data, as well as archival VLA 
data in the case of MJD 53511. 
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Fig. 6.— Top: dereddened GRS 1915+105 radio/MIR/X-ray SEDs, built with the 
MJD 53511, MJD 53280, and MJD 53851 RXTE /PCA+REXTE and IRS data (as 
well as archival VLA data in the case of MJD 53511), and fitted with the model 
phabs( diskir+gaussian +bbodyrad+bbodyrad+powerlaw). 

Bottom: dereddened 5.20 to 14.50 /im IRS spectra of GRS 1915+105, obtained on MJD 53511, 
MJD 53280, and MJD 53851, used to fit the radio/MIR/X-ray SEDs. Note that from the left to 
the right, both SEDs and spectra are displayed by decreasing levels of irradiation. 



